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Abstract 

This  study  investigates  the  effects  of  the  relative  humidity  (RH)  of  the  reactants  on  the  cell  performance  and  local  transport  phenomena  in  proton 
exchange  membrane  fuel  cells  with  parallel  and  interdigitated  flow  fields.  A  three-dimensional  model  was  developed  taking  into  account  the  effect 
of  the  liquid  water  formation  on  the  reactant  transport.  The  results  indicate  that  the  reactant  RH  and  the  flow  field  design  all  significantly  affect 
cell  performance.  For  the  same  operating  conditions  and  reactant  RH,  the  interdigitated  design  has  better  cell  performance  than  the  parallel  design. 
With  a  constant  anode  RH  =  100%,  for  lower  operating  voltages,  a  lower  cathode  RH  reduces  cathode  flooding  and  improves  cell  performance, 
while  for  higher  operating  voltages,  a  higher  cathode  RH  maintains  the  membrane  hydration  to  give  better  cell  performance.  With  a  constant 
cathode  RH=  100%,  for  lower  operating  voltages,  a  lower  anode  RH  not  only  provides  more  hydrogen  to  the  catalyst  layer  to  participate  in 
the  electrochemical  reaction,  but  also  increases  the  difference  in  the  water  concentrations  between  the  anode  and  cathode,  which  enhances  back- 
diffusion  of  water  from  the  cathode  to  the  anode,  thus  reducing  cathode  flooding  to  give  better  performance.  However,  for  higher  operating  voltages, 
the  cell  performance  is  not  dependent  on  the  anode  RH. 

©  2007  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Fuel  cells  have  attracted  much  attention  in  recent  years  as  a 
promising  alternative  to  traditional  internal  combustion  engines 
due  to  their  high  power  density  and  ultra-low  emissions.  Of  the 
various  types  of  fuel  cells,  the  proton  exchange  membrane  fuel 
cell  (PEMFCs)  operates  at  near-room  temperatures  and  is  con¬ 
sidered  to  be  a  good  choice  for  automotive  applications  [1-3]. 

Water  management  is  critical  to  the  PEMFC  performance.  It 
is  well  known  that  the  currently  used  commercial  polymer  elec¬ 
trolyte  membrane  (PEM)  in  PEMFCs  must  be  well  hydrated 
to  maintain  high-proton  conductivity,  which  is  proportional  to 
the  water  content  of  the  polymer  membrane.  Decreasing  water 
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content  in  the  polymer  membrane  leads  to  reduced  performance 
due  to  increases  of  the  voltage  loss  in  the  membrane.  To  avoid 
dehydration  of  the  polymer  electrolyte,  the  inlet  flows  are  usu¬ 
ally  injected  with  water  vapor.  In  addition,  the  electrochemical 
reactions  at  the  cathode  produce  water  vapor,  which  may  con¬ 
dense  and  accumulate  in  the  pores  of  the  porous  medium  and 
water  is  also  transported  from  the  anode  to  the  cathode  due  to 
electro-osmosis.  Excessive  liquid  water  can  cause  flooding  of 
the  pores  in  the  catalyst  layers  (CL)  and  the  gas  diffusion  lay¬ 
ers  (GDL),  thus  increasing  the  higher  mass  transfer  resistance 
to  the  reactant  flow.  Generally,  proper  flow  field  design  in  the 
bipolar  plates  is  an  effective  means  to  prevent  cathode  flooding 
and  increase  reactant  transport. 

Numerous  studies  have  presented  theoretical  models  to 
describe  the  water  transport  in  PEMFCs  [4—29].  The  earliest 
and  most  cited  models  were  developed  by  Springer  et  al.  [4] 
and  Bemardi  and  Vebrugge  [5,6].  Springer  et  al.  [4]  presented 
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Nomenclature 

a  chemical  activity  of  water  vapor 

A jT{fl  exchange  current  density  at  anode  (A  m-3) 

A jfc  exchange  current  density  at  cathode  (A  m-3 ) 

b  source  term  of  variable  <f> 

C  mass  fraction 

Cf  quadratic  drag  factor 

^porous  equivalent  surface  diameter  of  porous  media  (m) 
D  mass  diffusivity  (m2  s-1) 

F  Faraday  constant  (96,487  C  mol- 1 ) 

i  current  density  (A  m-2) 

I  average  current  density  in  the  fuel  cell  (A  m-2) 

Iy  local  current  density  in  the  y  direction  (A  m-2) 
j  current  density  (A  m-3) 

kc  water  vapor  condensation  rate  coefficient  (s- 1 ) 
ke  water  vapor  evaporation  rate  coefficient  (s- 1 ) 
kp  permeability  (m2) 

M  molecular  weight  (kg  mol- 1 ) 

P  pressure  (atm) 

P'  perturbed  pressure  in  a  control  volume 
R  universal  gas  constant  (8.314  Jmol-1  K-1) 
s  volume  ratio  occupied  by  liquid  water 

Sc  corrected  source  term  in  the  concentration  equa¬ 
tion 

Sj  source  term  in  the  phase  potential  equation 
.S'l  source  term  accounting  for  the  liquid  water  effect 
Su  source  term  in  the  x  momentum  equation 

Sv  source  term  in  the  y  momentum  equation 

Sw  source  term  in  the  z  momentum  equation 

S'  surface  area  (m2) 

t  time  (s) 

T  temperature  (K) 

u  x  direction  velocity  (ms-1) 

v  y  direction  velocity  (ms-1) 

Veen  operating  voltage  (V) 

V  volume  (m3) 

w  z  direction  velocity  (ms-1) 

x  x  direction  coordinate  (m) 

y  y  direction  coordinate  (m) 

Z  z  direction  coordinates  (m) 

Zf  species  valence 

Greek  letters 

aa  electrical  transfer  coefficient  in  the  forward  reac¬ 
tion 

ac  electrical  transfer  coefficient  in  the  backward 
reaction 
e  porosity 

0  dependent  variables 

0  phase  potential  function 

I I  overpotential  (V) 

X  membrane  water  content 

v  kinematic  viscosity  (m2  s- 1 ) 

p  density  (kg  m-3) 


crm  membrane  electrical  conductivity 
x  tortuosity  of  the  pores  in  the  porous  medium 

exchange  coefficient 

subscripts 
a  anode 

c  cathode 

channel  channel 
CL  catalyst  layer 

eff  effective 

GDL  gas  diffusion  layer 
H+  hydrogen  ion 

Hi  hydrogen 

H2O  water 

k  Mi  component  of  the  mixture 

MEM  membrane 

O2  oxygen 

porous  porous  medium 

sat  saturation 

total  total 

x  x  direction 

y  y  direction 

z  z  direction 

superscripts 

ref  reference  value 


a  one-dimensional  steady-state  model  for  a  PEMFC,  in  which 
water  diffusion  coefficients,  electro-osmotic  drag  coefficient, 
water  sorption  isotherms,  and  membrane  conductivity  were  mea¬ 
sured  as  a  function  of  the  membrane  water  content.  With  the 
assumption  of  a  fully  hydrated  membrane,  Bernardi  and  Ver- 
brugge  [5]  derived  a  mathematical  model  for  the  cathode  side 
of  a  PEMFC  and  found  that  the  water  transport  is  a  compli¬ 
cated  function  of  the  cell  operating  conditions.  They  further 
presented  a  full  cell  model  [6]  to  investigate  the  water  transport 
mechanisms. 

Other  important  early  fuel  cell  models  included  those  by 
Fuller  and  Newman  [7]  and  Nguyen  and  White  [8].  Fuller 
and  Newman  [7]  presented  a  quasi-two-dimensional,  along- 
the-channel  model  with  the  assumption  of  a  constant  water 
diffusivity  in  the  membrane  to  examine  the  water  management 
problems  in  a  PEMFC.  Nguyen  and  White  [8]  also  developed 
a  steady-state,  two-dimensional  model,  including  the  effect  of 
electro-osmosis,  water  diffusion,  and  heat  transfer  from  the  solid 
phase  to  the  gas  phase  and  the  latent  heat  due  to  water  evapora¬ 
tion  and  condensation,  to  investigate  the  effectiveness  of  various 
humidification  designs.  Okada  et  al.  [9]  studied  the  water  trans¬ 
port  on  the  anode  side  and  gave  a  linear  transport  equation  based 
on  the  water  diffusion  and  electro-osmotic  water  drag  to  analyze 
the  water  concentration  profiles.  Their  results  showed  that  the 
water  concentration  profile  is  influenced  by  the  membrane  thick¬ 
ness  and  humidity  and  the  current  density  in  the  cell.  Gurau  et 
al.  [10]  considered  the  variations  of  the  concentrations  and  par¬ 
tial  pressures  in  the  gas  channels  in  a  two-dimensional  model 
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for  the  entire  PEMFC  sandwich.  They  further  derived  a  half¬ 
cell  model  for  the  cathode  side  and  obtained  rigorous  analytical 
solutions  which  account  for  the  liquid  water  content  in  the  GDL 
[11]. 

Yi  and  Nguyen  [12]  developed  an  along-the-channel  model 
to  evaluate  the  effects  of  various  design  and  operating  parame¬ 
ters  on  the  performance  of  PEMFCs  to  show  that  humidification 
of  the  anode  reactants  is  required  to  enhance  the  membrane 
conductivity  and  liquid  injection  and  higher  humidification  tem¬ 
peratures  can  improve  cell  performance  by  introducing  more 
water  into  the  anode.  Hsing  and  Futerko  [13]  developed  a  two- 
dimensional  model  coupling  the  fluid  flow,  mass  transport  and 
electrochemistry  of  a  PEMFC  taking  into  account  the  varia¬ 
tion  of  the  liquid  water  diffusion  coefficient  in  the  membrane. 
Baschuk  and  Li  [14]  considered  the  effect  of  variable  degrees 
of  flooding  on  the  cathode  and  presented  a  model  to  evaluate 
its  influence  on  cell  performance.  Um  et  al.  [15]  presented  a 
transient  multidimensional  model  to  simulate  multi-component 
transport  in  a  PEMFC  and  explored  hydrogen  dilution  effects 
in  the  anode  feed.  Djilali  and  Lu  [16]  focused  on  the  model¬ 
ing  of  non-isothermal  and  non-isobaric  effects  and  performed 
an  analysis  of  the  cell  performance  and  water  transport  over 
a  range  of  operating  current  densities.  Kulikovsky  [17]  used  a 
quasi-three-dimensional  model  of  the  water  transport  in  PEM¬ 
FCs  to  investigate  the  non-linear  diffusion  of  liquid  water  in 
the  membrane.  Ge  and  Yi  [18]  developed  a  two-dimensional 
model  to  investigate  the  effects  of  operating  conditions  and 
membrane  thickness  on  the  water  transport.  They  used  an  effec¬ 
tive  porosity  to  account  for  the  effect  of  liquid  water  on  the  gas 
transport  to  simplify  the  two-phase  flow  model  in  the  porous  lay¬ 
ers.  Siegel  et  al.  [19]  derived  a  comprehensive  two-dimensional 
model  that  included  multi-component  and  multi-phase  transport 
both  along  the  gas  channel  and  through  the  membrane  elec¬ 
trode  assembly  (MEA).  Their  results  illustrated  the  importance 
of  water  transport  within  the  cell.  Mazumder  and  Cole  [20,21] 
developed  a  three-dimensional  model  to  predict  the  effect  of 
the  formation  of  liquid  water  on  the  performance  of  PEMFCs. 
Yan  and  co-workers  [22-32]  developed  two-dimensional  and 
three-dimensional  models  for  the  reactants  and  water  transport 
and  the  cell  performance  in  PEMFCs  with  various  flow  field 
designs.  The  effects  of  liquid  water  formation  on  the  reactant 
transport  were  taken  into  account  in  the  modeling  and  examined 
in  the  analysis.  Recently,  Wang  [33]  reviewed  the  fundamental 
models  for  PEMFCs. 

The  membrane  stability  is  determined  by  its  water  content. 
Proper  hydration  of  the  membrane  is  critical  for  maintain¬ 
ing  membrane  conductivity  and  mechanical  stability.  Various 
humidification  designs  such  as  internal  humidification,  external 
humidification  and  direct  injection  methods  are  used  in  PEM¬ 
FCs  to  maintain  hydration  of  the  polymer  membrane  [8,34 — 41], 
Recently,  Jang  et  al.  [42]  developed  a  two-dimensional  model 
to  investigate  the  effect  of  the  humidity  of  the  reactants  at  the 
inlet  on  the  detailed  gas  transport  and  cell  performance  of  the 
PEMFC  with  baffle-blocked  flow  field  designs.  Lee  and  Chu 
[43]  used  a  three-dimensional  model  to  investigate  the  loca¬ 
tion  of  the  gas-liquid  interface  for  various  cathode  humidity 
conditions  in  the  cathode  gas  diffusion  layer  for  the  conven¬ 


tional  flow  field.  In  their  study,  only  part  of  the  cell  was 
simulated  to  reduce  the  computational  effort.  These  studies 
revealed  that  the  relative  humidity  (RH)  critically  affects  cell 
performance. 

These  results  have  shown  that  the  optimal  cell  performance  is 
strongly  dependent  on  the  water  management.  Proper  humidifi¬ 
cation  of  the  reactants  can  improve  membrane  hydration  and  an 
appropriate  flow  field  design  in  the  bipolar  plates  can  increase 
reactant  transport  and  liquid  water  removal.  Though  some  stud¬ 
ies  have  analyzed  the  effect  of  the  relative  humidity  of  the 
reactants  on  the  PEMFC  performance,  most  previous  investiga¬ 
tions  have  utilized  simplified  two-dimensional  models  or  only 
considered  some  essential  parts  of  the  fuel  cell  in  the  compu¬ 
tational  domain  to  reduce  the  computational  time.  This  paper 
presents  a  three-dimensional,  full  cell  model  to  investigate  the 
effects  of  the  relative  humidity  of  the  reactants  on  both  the  anode 
and  cathode  sides  on  the  cell  performance  for  PEMFCs  with  par¬ 
allel  and  interdigitated  flow  fields.  The  effects  of  liquid  water 
formation  on  the  reactant  transport  are  taken  into  account  in 
the  model.  The  local  oxygen/hydrogen  mass  flow  rates,  liq¬ 
uid  water  distributions,  and  local  current  densities  in  the  cell 
for  both  flow  fields  are  examined  to  validate  the  simulation 
results. 

2.  Analysis 

A  three-dimensional  model  of  the  full  cell  was  developed  to 
analyze  the  electrochemical  reactions  and  transport  phenomena 
of  the  reactants  and  products  in  the  cell.  The  cell  was  divided  into 
the  anode  flow  channels,  membrane  electrode  assembly  (includ¬ 
ing  the  anode  GDL,  anode  CL,  proton  exchange  membrane, 
cathode  CL,  and  cathode  GDL),  and  cathode  flow  channels  as 
shown  in  Fig.  1.  The  governing  equations  include  the  mass, 
momentum,  species  and  electrical  potential  conservation  equa¬ 
tions.  The  model  assumes  that  the  system  is  three-dimensional 
and  steady;  the  inlet  fuels  are  ideal  gases;  the  system  is  isother¬ 
mal;  the  flow  is  laminar;  the  fluid  is  incompressible;  the  thermal 
properties  are  constant;  the  porous  layers  including  the  gas  dif¬ 
fusion  layer,  catalyst  layer  and  PEM,  are  isotropic.  The  present 
model  is  based  on  the  so-called  multiphase  multicomponent 
mixture,  originally  proposed  by  Wang  et  al.  The  transport  equa¬ 
tions  for  the  three-dimensional  PEM  fuel  cell  are: 


Fig.  1.  Schematic  of  the  cross-section  of  a  PEMFC  (not  to  scale  and  with  only 
two  channels  and  two  ribs  shown). 
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In  the  momentum  equations,  Su,  Sv  and  Sw  are  the  source 
terms  based  on  the  Darcy  drag  forces  in  the  x,  y,  and  z  direc¬ 
tions  imposed  by  the  pore  walls  on  the  fluid  which  usually 
cause  significant  pressure  drops  across  the  porous  media.  The 
details  of  Su,  Sv  and  Sw  for  the  various  layers  are  listed  in 
Table  1.  In  Table  1,  eeff  is  the  effective  porosity  of  the  porous 
material,  Cp  is  the  quadratic  drag  factor,  and  Zf  is  the  ion 
valence  in  the  PEM.  In  addition,  Dk,eS  is  the  effective  diffu¬ 
sion  coefficient  of  the  Ath  species,  which  was  calculated  using 
the  Bruggeman  equation  [44],  and  Ap  are  the  permeabilities 
of  the  porous  materials.  The  relationship  between  the  poros¬ 
ity  and  the  permeability  was  described  by  the  Blake-Kozeny 
equation  [45]: 


(5) 


where  dpoi ous  is  the  equivalent  pore  diameter  of  the  porous 
material  expressed  as  ^porous =  614 porous/ ■V/p0rnlls : 

•  Species  concentration  equation: 


(  dck  dck 


3CA 
W  dz) 


3  2Ck 
dy 2 


+  SC  +  SL  (6) 


where  Ck  is  the  concentration  of  the  Ath  component  and  Sc  and 
Sl  are  the  source  terms  for  the  chemical  reaction  and  the  liquid 
water  in  the  species  concentration  equation.  Sc  differs  for  the 
various  reactant  gases,  e.g.  Sc  is  — ya/2FCtotai,a  for  hydrogen, 
-jc/4FCtotal.c  for  oxygen,  and  jc/2FCtotai,c  for  water  vapor.  The 
parameters  ja  and  jc  denote  the  current  densities  on  the  anode  and 
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cathode  sides,  which  can  be  calculated  using  the  Butler- Volmer 
equations  [10]: 


(CA  1 

L(«a  F/RT)n  1 

[  e(acF/RT)r,\ 

(C?L \  1 

fetal  F/RT)n  1  1 

[<%)  1 

[  eCacF/sr^J 

(7) 

(8) 


where  A  jgef  is  the  reference  exchange  current  density,  aa  and  ac 
are  the  electrical  charge  transport  rates  in  the  anode  and  cathode 
catalyst  layers,  rj  the  overpotential,  F  the  Faraday’s  constant,  R 
the  ideal  gas  constant  and  T  is  the  fuel  cell  temperature.  The 
local  current  density  distribution  was  calculated  from  the  phase 
potential  equation: 


the  porous  media  cause  capillary  forces  to  dominate  the  liquid 
water  transport.  However,  the  actual  expression  of  this  force  can¬ 
not  be  formulated.  To  simplify  the  complex  two-phase  flow,  the 
simplified  liquid  water  transport  equation,  originally  proposed 
by  Mazumder  et  al.  [21],  was  used  in  the  present  model.  In  addi¬ 
tion,  the  liquid  water  formation  will  affect  the  effective  porosity 
and  the  mass  diffusivity  in  the  gas  diffusion  and  catalyst  lay¬ 
ers.  This  approach  was  also  adopted  by  Nguyen  et  al.  [44],  The 
source  term,  .S’l,  in  the  species  concentration  equation,  which  is 
a  function  of  the  liquid  water  saturation,  was  evaluated  as  [21] 

_  \  -MH,okcEeffC^0(PH2o  -  Fsat),  if  Ph20  >  Fsat 

Sl  =  <  pRi 

1  ke£ef{S(P sat  -  Fh2o)>  if  /JH20  <  Psat 

(17) 


where  is  the  phase  potential  and  erm  is  the  membrane  con¬ 
ductivity  which  is  strongly  dependent  on  the  membrane  water 
content  and  can  be  calculated  using  the  formulas  developed  by 
Springer  [4]: 


where  the  saturation  rate,  s,  was  defined  as  the  ratio  of  the  pore 
volume  occupied  by  liquid  water  to  the  total  pore  volume  in  the 
porous  media,  M  the  molecular  weight,  kc  and  ke  are  the  vapor 
condensation  and  evaporation  rate  coefficients,  and  Ps at  is  the 
water  vapor  saturation  pressure  calculated  using  [25,26,39]: 

psat  -  io_2  1794+0-02953:r_9-1837xlo_5:r2+1-4454x10”77'3  (18) 


».<r>  =  <exp[l268  (dj-i)] 


(10) 


The  effective  porosity  of  the  porous  media  was  then  modified 
to  account  for  the  liquid  water: 


where  <r^f  is  the  membrane  reference  conductivity  expressed  as 

<f  =  0.005139A.  -  0.00326  (11) 

0.043  +  17.81a  -  39.85a2  +  36.0a3,  0  <  a  <  1 

“  “  (12) 

14+  1.4(a-  1),  1  <  a  <  3 

where  a  is  the  water  activity  given  by  a  =  Cn2oRT/ Psat,  and  A. 
is  the  water  content  in  the  PEM.  Sj  in  Eq.  (9)  is  the  electrical 
source  term,  which  is  zero  in  the  PEM  without  electrochemical 
reactions  and  is  — ja  or  —jc  on  the  anode  or  cathode  sides.  The 
current  density,  i,  is  related  to  the  phase  potential,  0,  as 


fieff  =  e(l  -  s)  (19) 

Eqs.  (l)-(4),  (6)  and  (16)  form  a  complete  set  of  governing  equa¬ 
tions  for  the  PEMFC.  The  boundary  conditions  at  the  anode  flow 
channels  and  the  cathode  flow  channels  are  that  the  inlet  flow 
rates  are  constant,  the  inlet  gas  compositions  are  constant,  and 
the  flows  are  fully  developed  at  the  anode  and  cathode  flow  chan¬ 
nels  outlets.  The  solid  walls  are  no  slip  with  zero  flux  boundary 
conditions.  At  the  interfaces  between  the  gas  channels,  the  dif¬ 
fusion  layers,  the  catalyst  layers,  and  the  PEM,  the  velocities, 
mass  fractions,  momentum  fluxes,  and  mass  fluxes  are  assumed 
equal.  More  details  were  given  by  Yan  and  co-workers  [26]. 


30 


(13)  3.  Numerical  method 


30 


30 


(14) 

(15) 


Eq.  (9)  can  then  be  rewritten  as 

dix  3  iy  3  iz 

T  +  it  +  IT  =  SJ 

ox  ay  o  z 


(16) 


During  cell  operation,  the  partial  pressure  of  the  water  in  the 
electrode  may  exceed  its  saturation  pressure  if  the  local  water 
concentration  is  high.  Therefore,  when  the  water  vapor  partial 
pressure  was  greater  than  the  water  vapor  saturation  pressure, 
the  water  vapor  was  assumed  to  condense  and  fill  the  pores  in  the 
porous  media.  Operating  the  cell  at  high  reaction  rates  may  cause 
severe  mass  transport  overpotential  because  the  diffusing  com¬ 
ponents  are  blocked.  Furthermore,  the  extremely  small  pores  in 


The  SIMPLE  (semi  implicit  method  for  pressure-linked  equa¬ 
tion)  algorithm,  developed  by  Patankar  [46],  was  employed 
to  solve  the  governing  equations.  The  cell  performance 
was  simulated  for  a  miniature  cell  with  dimensions  of 
x  x  y  x  z  =  23  mm  x  23  mm  x  2.845  mm.  The  cell  had  12  chan¬ 
nels  (6  inlet  flow  channels  and  6  outlet  flow  channels  for  the 
interdigitated  flow  channel  design)  and  1 1  ribs,  all  1  mm  thick. 
The  diffusion  layer  was  0.4  mm  thick,  the  catalyst  layer  was 
0.005  mm  thick,  and  the  PEM  was  0.035  mm  thick  as  shown  in 
Fig.  2.  The  anode  flow  channels  were  assumed  to  be  in  parallel 
with  1  mm  channel  and  rib  widths  since  the  anode  flow  chan¬ 
nel  has  little  effect  on  the  cell  performance,  while  the  cathode 
flow  channels  were  assumed  to  be  parallel  or  interdigitated  flow 
channels  with  1  mm  channel  and  rib  widths.  The  base  operat¬ 
ing  conditions  for  the  fuel  cell  were  assumed  to  be  a  fuel  cell 
temperature  of  323  K,  the  reactant  gases  on  the  anode  side  were 
hydrogen  and  water  vapor  with  a  relative  humidity  of  100%,  the 
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flow  field. 

reactant  gases  on  the  cathode  side  contained  oxygen,  nitrogen, 
and  water  vapor  with  a  relative  humidity  of  100%,  the  inlet  flow 
rate  on  the  anode  side  was  260  cc  min-1,  the  inlet  flow  rate  on 
the  cathode  side  was  700  cc  min-1,  and  the  inlet  pressures  on 
the  anode  and  cathode  sides  were  both  1  atm.  All  of  the  fixed 
parameters  used  in  the  model  are  listed  in  Table  2.  To  study  the 
effect  of  the  relative  humidity  on  the  cathode  side,  RHC  was  cho¬ 
sen  as  25%,  50%,  75%,  and  100%  with  a  constant  RHa  =  100%. 


Table  2 

Fuel  cell  parameters 


Parameter 

Value 

A/fa 

9  x  108  Am“3 

A/f 

1.5  x  102  Am"3 

“a 

0.5 

1.5 

1 

^channel 

oo  m2 

«GDL 

0.5 

tgdl 

1.5 

&GDL 

1.76  x  10_1°  m2 

SCL 

0.4 

tcl 

1.5 

*CL 

1.76  x  10-11  m2 

SMem 

0.28 

TMem 

Dagan  Model 

^Mem 

1.8  x  10-18  m2 

/  (A  m  2) 


Fig.  3.  Influence  of  the  number  of  elements  on  the  polarization  curves. 


Similarly,  RHa  was  chosen  as  25%,  50%,  75%,  and  100%  with 
a  constant  RHC  =  100%  to  analyze  the  effect  of  the  RHa. 

The  grid  independence  was  examined  in  preliminary  test 
runs.  Four  non-uniformly  distributed  grid  configurations  were 
evaluated  for  the  parallel  flow  channel  design  with  the  num¬ 
bers  of  elements  in  the  x,  y  and  z  directions  being:  (I) 
70  x  70  x  25,  (II)  93  x  93  x  33,  (III)  116  x  116  x  33,  and  (IV) 
116x  116x41. The  influence  of  the  number  of  elements  on  the 
polarization  curves  (/-Vceii  curves)  for  the  fuel  cell  is  shown  in 
Fig.  3.  Grid  (III)  was  chosen  for  the  simulations  as  a  tradeoff 
between  accuracy  and  execution  time.  The  coupled  set  of  equa¬ 
tions  was  solved  iteratively  and  the  solution  was  considered  to 
be  converged  when  the  relative  error  in  each  field  between  two 
consecutive  iterations  was  less  than  10-6. 

The  numerical  results  were  validated  by  comparing  the 
present  predictions  with  previous  experimental  results.  Fig.  4 
compares  the  calculated  polarization  curve  with  experimental 
data  [47]  for  a  fuel  cell  with  the  parallel  flow  channel  design 
and  an  area  of  140  mm  x  140  mm.  There  is  only  a  small  differ¬ 
ence  between  the  numerical  results  and  the  experimental  data. 
Hence,  the  three-dimensional  numerical  model  can  be  used  to 
accurately  analyze  the  effects  of  the  flow  channel  area  ratio  and 
cathode  fuel  flow  rate  on  the  cell  performance. 


/  (A  m  2) 


Fig.  4.  Comparison  of  predicted  and  measured  polarization  curves  for  a  fuel  cell 
with  the  parallel  flow  channel  design. 
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4.  Results  and  discussion 

The  PEMFC  performance  is  strongly  dependent  on  the  water 
management  in  the  cell  which  includes  two  aspects.  First,  the 
membrane  hydration  must  be  maintained  for  proton  transport, 
since  the  ohmic  resistance  of  the  membrane  increases  as  the 
membrane  water  content  decreases,  which  reduces  cell  perfor¬ 
mance.  Secondly,  the  electrochemical  reaction  on  the  cathode 
side  produces  water  vapor  and  when  the  partial  pressure  of 
the  water  vapor  is  higher  than  its  saturation  pressure,  the  water 
vapor  condenses  to  form  liquid  water.  In  addition,  water  is  also 
transported  from  the  anode  to  the  cathode  by  electro-osmotic 
drag.  When  the  water  generation  rate  on  the  cathode  side  by 
electro-osmotic  drag  and  the  oxygen  reduction  reaction  exceeds 
the  water  removal  rate  from  the  cathode  by  back-diffusion  to 
the  anode,  evaporation,  and  capillary  transport  of  liquid  water 
through  the  cathode  gas  diffusion  layer  and  catalyst  layer,  the 
cathode  becomes  flooded.  Excessive  accumulated  liquid  water 
blocks  the  gas  pores  in  the  diffusion  layer  and  the  catalyst  layer, 
which  forms  a  barrier  over  the  catalyst  active  surface  in  the 
catalyst  layer,  worsening  the  performance.  In  general,  humid¬ 
ification  of  the  reactants  can  ensure  membrane  hydration,  but 
may  induce  cathode  flooding.  Proper  flow  field  design  is  con¬ 
sidered  to  be  an  effective  means  to  prevent  cathode  flooding 
because  it  can  increase  reactant  transport  rates  and  enhance  liq¬ 
uid  water  removal.  This  discussion  analyzes  the  effects  of  the 
relative  humidity  of  the  reactants  on  both  the  cathode  and  anode 
sides  on  the  cell  performance  for  parallel  and  interdigitated  flow 
fields. 

4.1.  Effect  of  cathode  relative  humidity 


I  (A  m  2) 


Fig.  5.  Effect  of  cathode  relative  humidity  on  the  polarization  curves:  (a)  parallel 
flow  field;  (b)  interdigitated  flow  field. 


Fig.  5(a)  and  (b)  show  the  I-V ceii  polarization  curves  for  par¬ 
allel  and  interdigitated  flow  fields  for  various  cathode  relative 
humidities  with  a  constant  RHa  =  100%.  At  low-operating  volt¬ 
ages  for  both  flow  fields,  the  cell  performance  decreases  as  RHC 
increases.  The  cell  performance  is  the  best  for  RHC  =25%  and 
the  worst  for  RHC  =  100%  because  at  low-operating  voltages,  the 
electrochemical  reactions  are  stronger  with  more  liquid  water 
produced  on  the  cathode  side,  so  the  back-diffusion  of  water 
and  the  high  anode  relative  humidity,  RHa  =  100%,  provide  suf¬ 
ficient  water  for  the  membrane.  Thus,  the  cell  performance  is 
mainly  dependent  on  the  cathode  mass  transport  limitations  due 
to  the  liquid  water  blockage  effect.  When  RHC  is  lower,  the 
oxygen  concentration  in  the  reactants  is  higher  and  the  water 
vapor  concentration  on  the  cathode  side  is  lower,  which  reduces 
the  cathode  flooding,  so  the  cell  performance  increases.  How¬ 
ever,  as  the  operating  voltage  increases,  the  effect  of  RHC  on  the 
cell  performance  reverses  (this  operating  voltage  is  referred  to 
here  as  the  reversal  operating  voltage).  The  cell  performance  is 
the  worst  for  RHC  =25%,  better  for  RHC  =  50%,  and  the  best 
for  RHC  =  75%  with  little  difference  for  RHC  =  100%.  These 
changes  occur  because  as  the  operating  voltage  increases,  the 
electrochemical  reactions  slow  and  only  a  small  amount  of  liq¬ 
uid  water  is  produced  on  the  cathode  side,  so  cathode  flooding  is 
no  longer  the  key  factor  affecting  cell  performance.  As  the  water 
concentration  on  the  cathode  side  decreases,  the  mass  flow  rate  of 


water  transport  from  the  cathode  to  the  anode  by  back-diffusion 
decreases;  thus,  the  membrane  hydration  gradually  becomes  the 
key  factor  affecting  cell  performance.  Higher  RHC  increase  the 
water  concentration  on  the  cathode  side  and  the  back-diffusion 
of  water,  which  provides  more  water  for  the  membrane,  result¬ 
ing  in  better  cell  performance.  Comparison  of  Fig.  5(a)  and  (b) 
indicates  that  at  the  same  operating  conditions,  the  flow  field 
design  also  significantly  affects  the  cell  performance.  For  operat¬ 
ing  voltages  lower  than  0.7  V,  the  interdigitated  flow  field  design 
has  better  cell  performance  than  the  parallel  flow  field  design. 
Furthermore,  the  reversal  operating  voltages  for  the  parallel  and 
interdigitated  flow  field  designs  also  differ.  For  the  parallel  flow 
field  design,  the  reversal  operating  voltage  is  about  0.7  V,  while 
for  the  interdigitated  flow  field  design  it  is  about  0.575  V.  The 
baffles  at  the  end  of  the  inlet  flow  channels  in  the  interdigitated 
flow  field  change  the  reactant  transport  into  the  GDF  and  CL 
from  diffusion  to  forced  convection.  This  forced  convection  not 
only  leads  to  more  reactant  transport  into  the  GDL  and  CL  for 
the  electrochemical  reaction,  but  also  induces  larger  shear  forces 
which  increase  liquid  water  removal.  Therefore,  the  interdigi¬ 
tated  flow  field  design  has  better  performance  than  the  parallel 
flow  field  design  at  low-operating  voltages.  Additionally,  cath¬ 
ode  flooding  in  the  interdigitated  flow  field  design  is  significantly 
reduced  compared  with  the  parallel  flow  field  design,  so  that  the 
mass  flow  rate  of  water  from  the  cathode  to  the  anode  by  back- 
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middle  cross-section  in  the  PEM  at  operating  voltages  of  0.3  and  0.7  V:  (a) 
parallel  flow  field;  (b)  interdigitated  flow  field. 


the  cathode  GDL-CL  interface  at  operating  voltages  of  0.3  and  0.7  V:  (a)  parallel 
flow  field;  (b)  interdigitated  flow  field. 


diffusion  decreases  and  the  membrane  hydration  has  a  more 
significant  effect  on  cell  performance  at  lower  operating  volt¬ 
ages  than  with  the  parallel  flow  field  design.  Thus,  the  reversal 
operating  voltage  for  the  interdigitated  flow  field  design  is  lower 
than  for  the  parallel  flow  field  design. 

Fig.  6(a)  and  (b)  show  the  effect  of  the  cathode  relative 
humidities  on  the  local  current  densities  on  the  middle  cross- 
section  (at  y  =  1 1 .5  mm)  in  the  PEM  at  operating  voltages  of  0.3 
and  0.7  V  for  the  parallel  and  interdigitated  flow  field  designs 
with  a  constant  RHa  =  100%.  For  the  parallel  flow  field  design, 
the  maximum  current  densities  occur  at  the  center  of  the  flow 
channels  with  the  minimum  current  densities  at  the  center  of 
the  ribs,  while  for  the  interdigitated  flow  field  design,  the  maxi¬ 
mum  current  densities  occur  at  the  inlet  flow  channels  with  the 
minimum  current  densities  at  the  outlet  flow  channels.  At  the 
operating  voltage  of  0.3  V,  for  both  flow  field  designs,  the  min¬ 
imum  and  maximum  local  current  densities  for  RHC  =25%  are 
both  higher  than  for  the  other  three  cathode  relative  humidities, 
so  the  cell  performance  is  best  for  RHC  =  25%.  However,  at  the 
operating  voltage  of  0.7  V,  for  both  flow  field  designs,  the  local 
current  densities  for  RHC  =  25  %  are  lower  than  for  the  other  three 
cathode  relative  humidities,  so  the  cell  performance  is  worst  for 
RHC  =  25%. 

Fig.  7(a)  and  (b)  show  the  effect  of  the  cathode  relative 
humidities  on  the  oxygen  mass  flow  rates  across  the  cathode 


GDL-CL  interface  (y=  11.5  mm)  at  operating  voltages  of  0.3 
and  0.7  V  for  the  parallel  and  interdigitated  flow  field  designs 
with  a  constant  RHa  =  100%.  The  oxygen  mass  flow  rate  across 
the  GDL-CL  interface  represents  the  amount  of  oxygen  entering 
the  CL  per  unit  time;  thus,  higher  oxygen  mass  flow  rates  mean 
that  more  oxygen  participates  in  the  electrochemical  reaction  in 
the  CL,  resulting  in  larger  local  current  densities.  Comparison  of 
Figs.  6  and  7  indicates  that  the  local  current  densities  correlate 
perfectly  with  the  oxygen  mass  flow  rates,  which  validates  the 
present  model. 

Fig.  8(a)  and  (b)  show  the  effect  of  the  cathode  relative 
humidities  on  the  liquid  water  concentrations  along  the  cathode 
GDL-CL  interface  (y=  11.5  mm)  at  operating  voltages  of  0.3 
and  0.7  V  for  the  parallel  and  interdigitated  flow  field  designs 
with  a  constant  RHa=  100%.  For  both  flow  field  designs,  the 
liquid  water  concentrations  at  the  lower  operating  voltage  of 
0.3  V  are  far  higher  than  at  the  higher  operating  voltage  of 
0.7  V,  which  supports  the  observation  that  the  mass  transport 
limitation  due  to  cathode  flooding  is  the  main  factor  affecting 
the  cell  performance  at  lower  operating  voltages.  However,  at 
higher  operating  voltages,  the  liquid  water  concentration  at  the 
cathode  is  lower,  so  concentration  polarization  has  little  effect 
on  the  cell  performance.  Although  the  anode  relative  humidity 
is  100%,  the  mass  flow  rate  of  water  from  the  cathode  to  the 
anode  by  back-diffusion  decreases  at  higher  operating  voltages; 
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along  the  cathode  GDL-CL  interface  at  operating  voltages  of  0.3  and  0.7  V:  (a) 
parallel  flow  field;  (b)  interdigitated  flow  field. 
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Fig.  9.  Effect  of  anode  relative  humidity  on  the  polarization  curves:  (a)  parallel 
flow  field;  (b)  interdigitated  flow  field. 


thus,  the  membrane  hydration  gradually  becomes  the  key  factor 
affecting  the  cell  performance.  Fig.  8  shows  that  at  the  lower 
operating  voltage  of  0.3  V,  for  both  flow  field  designs,  the  liq¬ 
uid  water  concentrations  are  largest  for  RHC  =  100%  with  the 
severest  concentration  polarization,  so  the  cell  performance  is 
the  worst.  Although  at  the  higher  operating  voltage  of  0.7  V, 
for  both  flow  field  designs,  the  liquid  water  concentrations  are 
also  largest  for  RHC  =  100%,  the  concentration  polarization  is 
no  longer  the  main  factor  affecting  the  cell  performance.  Higher 
liquid  water  concentrations  enhance  the  back-diffusion  of  water 
and  effectively  ensure  membrane  hydration,  which  improves 
cell  performance.  Comparison  of  Figs.  7  and  8  shows  that  the 
liquid  water  distributions  are  exactly  reversed  from  the  local 
oxygen  mass  flow  rate  distributions  at  the  lower  operating  volt¬ 
age  of  0.3  V.  For  the  parallel  flow  field  design,  the  maximum 
liquid  water  concentrations  occur  at  the  center  of  the  ribs  with 
the  minimums  at  the  center  of  the  flow  channels;  while  for  the 
interdigitated  flow  held  design,  the  maximum  liquid  water  con¬ 
centrations  occur  at  the  outlet  flow  channels  with  the  minimums 
at  the  inlet  how  channels.  This  occurs  because  as  the  liquid 
water  concentrations  increase,  more  pores  are  blocked  in  the 
GDL  and  CL  which  increases  the  oxygen  transport  resistance, 
resulting  in  decreased  oxygen  mass  how  rates.  Comparison  of 
Fig.  8(a)  and  (b)  indicates  that  the  liquid  water  concentrations 
for  the  interdigitated  how  held  design  are  far  less  than  for  the 


parallel  How  held  design,  which  conhrms  that  the  interdigitated 
how  held  design  has  more  efficient  liquid  water  removal  than 
the  parallel  how  held  design  and  signihcantly  reduces  cathode 
hooding.  Therefore,  as  shown  in  Fig.  5  the  concentration  polar¬ 
ization  in  the  polarization  curves  for  the  interdigitated  how  held 
design  is  not  significant  compared  to  that  in  the  paraliel  how  held 
design. 

4.2.  Effect  of  anode  relative  humidity 

Fig.  9(a)  and  (b)  show  the  /-Vceii  polarization  curves  for 
the  parallel  and  interdigitated  how  helds  at  various  anode  rel¬ 
ative  humidities  with  a  constant  RHC  =  100%.  For  both  how 
held  designs,  at  operating  voltages  greater  than  0.7  V,  the  /-Fceii 
polarization  curves  for  the  various  RHa  almost  coincide,  indicat¬ 
ing  that  the  cell  performance  is  not  dependent  on  RHa.  However, 
at  operating  voltages  lower  than  0.7  V,  the  cell  performance 
curves  for  the  various  RHa  start  to  differ.  As  the  RHa  decreases, 
the  limiting  current  densities  increase  and  the  cell  performance 
improves  because  the  molecular  weight  of  hydrogen  is  far  less 
than  that  of  water,  so  as  the  RHa  decreases  the  hydrogen  con¬ 
centration  in  the  reactants  is  signihcantly  increased.  In  addition, 
at  lower  operating  voltages,  the  reduced  RHa  increases  the  dif¬ 
ferences  in  the  water  concentrations  between  the  cathode  and 
anode,  which  not  only  reduces  cathode  hooding,  but  also  pro- 
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middle  cross-section  in  the  PEM  at  operating  voltages  of  0.3  and  0.7  V:  (a) 
parallel  flow  field;  (b)  interdigitated  flow  field. 

vides  more  water  for  the  membrane,  so  the  cell  performance 
improves. 

Fig.  10(a)  and  (b)  show  the  effect  of  anode  relative  humidi¬ 
ties  on  the  local  current  densities  on  the  middle  cross-section 
(at  y  =  1 1 .5  mm)  in  the  PEM  at  the  operating  voltages  of  0.3  and 
0.7  V  for  the  parallel  and  interdigitated  flow  field  designs  with 
a  constant  RHC  =  100%.  At  the  operating  voltage  of  0.7  V,  for 
both  flow  field  designs,  the  curves  of  the  local  current  density 
for  various  RHa  almost  coincide,  indicating  that  the  cell  per¬ 
formance  is  not  dependent  on  RHa.  At  the  operating  voltage 
of  0.3  V,  the  local  current  density  distributions  for  the  vari¬ 
ous  RHa  are  similar  to  those  in  Fig.  6  for  the  various  RHC. 
For  both  flow  held  designs,  the  local  current  densities  reach 
maximums  at  RHa  =  25%,  so  the  cell  performance  there  is 
best. 

Fig.  1 1  (a)  and  (b)  show  the  effect  of  the  anode  relative  humidi¬ 
ties  on  the  hydrogen  mass  flow  rates  across  the  anode  GDL-CL 
interface  (y=  11.5  mm)  at  operating  voltages  of  0.3  and  0.7  V 
for  the  parallel  and  interdigitated  flow  held  designs  with  a  con¬ 
stant  RHC  =  100%.  Because  a  change  in  the  RHa  signihcantly 
affects  the  hydrogen  concentration  in  the  reactants,  as  RHa 
decreases  the  hydrogen  mass  How  rates  signihcantly  increases 
since  there  is  more  hydrogen  transport  into  the  CL  to  participate 
in  the  electrochemical  reactions  which  improves  cell  perfor¬ 
mance.  Fig.  12(a)  and  (b)  show  the  effect  of  the  anode  relative 
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Fig.  1 1 .  Effect  of  anode  relative  humidity  on  the  hydrogen  mass  flow  rates  across 
the  anode  GDL-CL  interface  at  operating  voltages  of  0.3  and  0.7  V:  (a)  parallel 
flow  field;  (b)  interdigitated  flow  field. 


humidities  on  the  liquid  water  concentrations  along  the  cathode 
GDL-CL  interface  (y=  11.5  mm)  at  operating  voltages  of  0.3 
and  0.7  V  for  the  parallel  and  interdigitated  flow  held  designs 
with  a  constant  RHC  =  100%.  For  both  flow  held  designs,  the  liq¬ 
uid  water  concentrations  are  highest  at  RHa  =  100%,  and  as  RHa 
decreases  the  liquid  water  concentrations  gradually  decrease, 
indicating  that  back-diffusion  of  liquid  water  from  the  cathode 
to  the  anode  increases  as  the  anode  relative  humidity  decreases 
because  the  decreasing  RHa  increases  the  differences  between 
the  water  concentrations  at  the  cathode  and  anode.  Note  that 
previous  investigations  [34,35,41]  indicated  that  lower  anode 
relative  humidities  usually  resulted  in  a  higher  membrane  resis¬ 
tance  and,  thus;  higher  ohmic  voltage  losses  in  the  membrane 
because  the  cathode  relative  humidities  were  maintained  at  val¬ 
ues  far  less  than  100%  in  those  investigations.  However,  in  the 
present  paper,  the  cathode  humidity  is  maintained  at  100%;  thus, 
the  high  cathode  humidity  provides  sufficient  water  to  maintain 
membrane  hydration  by  back-diffusion  of  water.  Particularly, 
at  lower  operating  voltages,  the  concentration  polarization  is 
the  main  factor  affecting  the  cell  performance,  with  high-liquid 
water  concentrations  on  the  cathode  side  increasing  the  oxygen 
transport  resistance.  Because  a  decreasing  RHa  reduces  cathode 
flooding,  which  allows  more  oxygen  to  enter  the  CL  to  partici¬ 
pate  in  the  electrochemical  reaction;  thus,  the  cell  performance 
improves. 
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Fig.  12.  Effect  of  anode  relative  humidity  on  the  liquid  water  concentrations 
along  the  cathode  GDL-CL  interface  at  operating  voltages  of  0.3  and  0.7  V:  (a) 
parallel  flow  field;  (b)  interdigitated  flow  field. 

In  addition,  comparison  of  Figs.  5  and  9  shows  that  for  the 
same  inlet  mass  flow  rates  of  reactants  on  the  cathode  and  anode 
sides,  the  reduction  of  RHa  for  a  constant  RHC  =  100%  is  a  more 
effective  means  for  improving  cell  performance  than  reducing 
RHC  with  a  constant  RHa=  100%.  First,  both  RHa  =  100%  or 
RHC  =  100%  provide  good  membrane  hydration.  Secondly,  at 
lower  operating  voltages,  a  large  amount  of  liquid  water  is 
produced  by  the  cathode  electrochemical  reaction  with  the  back- 
diffusion  of  water  further  ensuring  membrane  hydration.  Thus, 
cathode  flooding  is  the  main  factor  affecting  cell  performance 
at  lower  operating  voltages.  The  reduction  of  RHa  can  lead  to 
more  water  transport  from  the  cathode  to  the  anode  by  back- 
diffusion,  which  reduces  cathode  flooding  and  mass  transport 
losses.  Decreasing  RHa  also  significantly  increases  the  hydro¬ 
gen  concentration  in  the  reactants  due  to  the  small  molecular 
weight  of  hydrogen. 

5.  Conclusions 

The  proper  relative  humidity  of  the  reactants  and  the  proper 
flow  field  designs  in  the  bipolar  plates  both  effectively  improve 
water  management  and  cell  performance  in  PEMFC.  This  paper 
presented  a  three-dimensional  numerical  model  based  on  the 
finite  volume  method  to  analyze  the  effects  of  the  relative 
humidities  of  the  reactants  on  both  the  cathode  and  anode  sides 


on  the  cell  performance  for  parallel  and  interdigitated  flow  fields . 

The  conclusions  drawn  from  the  analyses  are: 

(1)  At  the  same  operating  conditions  and  reactant  relative 
humidities,  the  interdigitated  flow  field  design  has  better 
cell  performance  than  the  parallel  flow  field  design  because 
the  baffle  changes  the  reactant  transport  into  the  GDL  and 
CL  from  diffusion  to  forced  convection,  which  increases 
reactant  transport  rates  and  enhances  liquid  water  removal. 

(2)  With  a  constant  RHa  =  100%  for  lower  operating  voltages, 
a  lower  RHC  reduces  cathode  flooding,  which  improves  cell 
performance.  For  higher  operating  voltages,  a  higher  RHC 
improves  membrane  hydration  which  gives  better  cell  per¬ 
formance. 

(3)  With  a  constant  RHC  =  100%  for  lower  operating  voltages, 
a  lower  RHa  not  only  provides  more  hydrogen  into  the 
CL  to  participate  in  the  electrochemical  reaction,  but  also 
increases  the  differences  between  the  water  concentrations 
at  the  anode  and  cathode,  which  enhances  the  back-diffusion 
of  water  from  the  cathode  to  the  anode;  thus,  reducing  the 
cathode  flooding.  As  a  result,  a  lower  RHa  gives  better  per¬ 
formance  than  a  higher  RHa.  For  higher  operating  voltages, 
the  cell  performance  is  not  dependent  on  RHa. 

(4)  For  the  same  reactant  inlet  flow  rates  on  the  cathode 
and  anode  sides,  the  reduction  of  RHa  with  a  constant 
RHC  =  100%  is  a  more  effective  means  to  improve  cell  per¬ 
formance  than  reducing  RHC  with  a  constant  RHa=  100%. 
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